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A topological insulator is characterized by a dichotomy between the interior and the edge of a 
finite system: While the bulk has a non-zero energy gap, the edges are forced to sustain excitations 
traversing these gaps. Originally proposed for electrons goverened by quantum mechanics, it has 
remained an important open question if the same physics can be observed for systems obeying 
Newton’s equations of motion. Here, we report on measurements that characterize the collective 
behavior of mechanical oscillators exhibiting the phenomenology of the quantum spin hall effect. The 
phononic edge modes are shown to be helical and we demonstrate their topological protection via 
the stability against imperfections. Our results open the door to the design of topological acoustic 
meta-materials that can capitalize on the stability of the surfaces phonons as reliable wave guides. 


The experimental hallmarks of the quantum ^in Hall 
effect (QSHE) in semi-conductor quantum wells^are two 
counter-propagating edge modes that differ by their spin 
degree of freedom. As long as time reversal symmetry 
is preserved, these two modes are independent and do not 
scatter into each other. Much of the interest in condensed- 
matter research involving topological states is driven by 
the use of these protected edge modes for technological 
applications such as spintronicsj^ magnetic device^ or 
quantum information processing.121 The transfer of the phe¬ 
nomenology of the QSHE from the quantum mechanical 
realm to classical mechanical systems is therefore both of 
fundamental interest and a gateway towards new design 
principles in mechanical meta-materials. 

Several key problems in the engineering of acoustic ma¬ 
terials can potentially be addressed by capitalizing on 
the physics of the QSHE. The edge channels are ro¬ 
bust c ounte r-parts to the well known whispering gallery 
modes.liSIIIl Any application that requires energy to be 
confined to the surfaces of some device, for example vi¬ 
bration insulators, can potentially make use of such edge 
states. In contrast to the whispering gallery modes, which 
are extremely sensitive to the shape of the surface,^ the 
topological edge modes are stable under a broad range 
of perturbations. Moreover, due to the stability of these 
modes, scattering-free phonon wave-guides of almost ar¬ 
bitrary shape can be realized. This in turn enables the 
engineering of robust acoustic delay lineJi^®^ useful for 
acoustic leasing 

How a mechanical system described by Newton’s equa¬ 
tions can reproduce the phenomenology of a quantum me¬ 
chanical model such as the QSHE has, however, remained 
an open question.E^ti^ Moreover, with a view to potential 
future applications, it is essential to demonstrate that the 
physics of the non-interacting and damping-free QSHE can 
be observed in a real mechanical system, which necessarily 
suffers from losses, disorder and non-linearities inherent to 
coupled mechanical oscillators. 

A quantum-mechanical lattice problem is described by 
a Schrodinger equation of the form 

ih-ipf = ( 1 ) 

where are the wave-function amplitudes for an electron 


with spin a on lattice site i and is the Hamiltonian 
matrix describing the QSHE.I^^HU Qn the other hand, the 
dynamics of a collection of classical harmonic oscillators is 
described by Newton’s equation of motion 

Xi — T^ijXj, ( 2 ) 

where Xi are the coordinates of N pendula {i = I,..., N) 
and Vij is the dynamical matrix containing the couplings 
between them. For either system, the existence and prop¬ 
erties of the edge modes are features of the eigenstates of 
PL OT V alone and depend neither on the interpretation of 
ipf versus Xi, nor on the nature of the dynamics, i.e., idt 
versus d^. Hence, we aim to design a dynamical matrix D 
incorporating the properties of the QSHE. 


a b c 



FIG. 1: Setup, a, Illustration of two pendula making up one 
effective site of our lattice model. Six springs are attached to 
each pendulum for the coupling to its nearest neighbors, b, 
Schematic (top) view of the different types of couplings. Each 
green rectangle depicts a logical site hosting an “a;” and a “y” 
pendulum swinging in the direction of the arrows. The longi¬ 
tudinal couplings are simple springs connecting the pendula. 
Transverse couplings are achieved via one or two lever arms, 
depending on the required sign of the coupling (see App. 0.C, 
Our 9 X 15 lattice made of 270 pendula. The two helical edge 
states known from the quantum spin Hall effect are indicated 
in blue and red. 
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FIG. 2: Helical edge mode dispersion, a, Time traces of two pendula. These can be interpreted as a two-dimensional 
trajectory, as shown. Steady states are displayed by colored disks representing their polarization. The radius of the circle 
corresponds to the mean deflection A and the black line indicates the position of the pendula at a given fixed time, b & c, 
Measured steady states at two different frequencies (2.114 and 2.380Hz). The circles are normalized to the strongest deflection. 
The excluded site at the centre of the bottom row is excited with left-circular polarization, d, Mean response of the system 
(average A) as a function of the excitation frequency showing an overall bandwidth of ~ 1.2 Hz. Error bars (smaller than the 
symbol size) are explained App. e, The mean response of the edge relative to the bulk. In the frequency ranges shaded in 
grey, the bulk response dominates. The frequencies marked with dashed lines correspond to panels b &: c. f, For edge-dominated 
modes, the evolution of the angle indicated by the black line dehnes a wave vector k. g, The frequency of the edge modes as 
a function of the wave vector. The bulk bands are shaded in gray. The color labels the polarization as before, establishing the 
helical nature of the edge excitations. The shape of the dark dots indicates their error (App. 0 while the red (blue) shaded region 
marks a 0.04 Hz band corresponding to the loss-induced broadening of the eigen-frequencies. 


We start from two independent copies of the Hofstadter 
modeP^ on the square lattice with flux $ = per pla- 
quette 

H =f ^ \r, s, a) (r, s ± 1, a| -b |r, s, a) {r ± 1, s, . 

r,s,a=± 

( 3 ) 

Here, (r, s) denotes the location on the lattice of size Lr x 
Ls, 4’s = 1‘1’ls and / represents the hopping amplitude. 
The two copies are labelled by a pseudo-spin index a = ±. 


H is symmetric under time reversal T and has three doubly 
degenerate bands, which are separated by non-zero gaps 
with a non-trivial topological Z 2 index, cf. App 
Therefore, in a finite system we expect one chiral edge 
state per pseudo-spin in both gaps. 

To take the step from the Hamiltonian "H to the dynam¬ 
ical matrix T), we need to ensure two properties. First, the 
eigenvalues of V are the squares of the eigenfrequencies of 
the dynamical system, cf. Eq. <§• Therefore, V has to 
be positive definite, which can always be achieved by the 
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FIG. 3: Beam splitter, a, Geometry of the beam splitter. 
An arbitrary polarization is injected at the indicated site at 
2.123 Hz. The relative weights A on the two shaded boxes 
(R/L) define a splitting ratio, b, Splitting ratio on the Bloch 
sphere. The north (south) pole of the sphere corresponds to left 
(right) circular polarizations of the drive and the color indicates 
the splitting ratio. The colormap is an interpolation based 
on 7000 data points. The maximal splitting is not reached 
at the poles, indicating the presence of disorder. The dots 
mark measurements along a great circle through the points of 
maximal splitting, c, Splitting ratio along the great circle. The 
black line marks a cosine expected for an optimal beam-splitter. 
The maximally reached splitting is 99.80(4)%. 


addition of a diagonal term oc /i. Second, V has to be 
strictly non-imaginary, as it encodes ideal springs. Owing 
to the time-reversal symmetry, this can be attained by a 
unitary transformation combining the local Kramers pairs 
a = ± 



U — U 0 1 lattice- 


these we obtain the mean deflection Ar^s as well as the 
polarization (left handed, linear or right handed) related 
to the lag between the x and y pendulum (see App. 0, 
cf. Fig.[|;a). 

To establish the existence of the edge states, we scan the 
frequency at which one edge site is excited. For every fre¬ 
quency we wait until the system reaches a steady state be¬ 
fore we extract the total response x = X)r s ^r,s/N, where 
N is the number of sites. Figure [^d) shows that the sys¬ 
tem responds appreciably between 1.7Hz and 2.9 Hz, as 
expected (see App. 0. Moreover, there are two regions 
with sequences of pronounced peaks. The width of these 
peaks is F 0.04 Hz, indicating the damping of the oscil¬ 
lators. In order to make the connection to edge modes we 
separate the system into the outermost line of lattice sites 
Xe = Eedge and the rest Xb = Ebuik Ar,s/A^buik- 

The relative weight XeliXb + Xe) is shown in Fig. [^e). 
There are three bands where the response lies mainly in 
the bulk which are separated by two frequency regions 
where the response is dominated by the edge. To illus¬ 
trate this further, we show in Fig. [^b) & (c) the recorded 
mode structure at a bulk and an edge frequency, respec¬ 
tively. 

Given the analogy to the QSHE effect, we expect the 
edge states to be helical. The edge spectrum uj{k), where 
k is the wave vector along the edge, can be extracted from 
the steady states shown in Fig. Beyond A^^s and 

the polarization, we determine the position of each two- 
dimensional oscillator at a given time [Xr,sito),yr,sito)]. 
The angle of this vector with respect to the positive x- 
direction defines a local phase cf. Fig. I^a). The 
change 


k — 4>r,s+l 4’r,s 


Each of the new coordinates Xr,s and yr,.s describes a one¬ 
dimensional oscillator at lattice site (r, s). In the follow¬ 
ing, we interpret {xr,s,yr,s) as the coordinates of a two- 
dimensional pendulum. It is then evident from the struc¬ 
ture of u that the eigenmodes of the system characterized 
by a = ± correspond to left and right circularly polarized 
motions. Hence, these two polarizations replace the notion 
of pseudo-spin up and down of the quantum mechanical 
problem. 

We implement the dynamical matrix D = WT-LU using 
pendula as depicted in Fig. V inherits its structure 
from the Landau gauge in Eq. In s-direction, there 
is no cross-coupling between Xr,s and 2/r,s±i (no Peierls 
phase in "H). Along the r-direction, T) also mixes the x- 
and j/-degrees of freedom with a strength that depends on 
the r-coordinate (see App. 0. The total system consists 
of 270 pendula in a lattice of Lr x Lg = 9 x 15 sites. 
The bare eigenfrequencies of the pendula due to gravity is 
a;o/27r « 0.75 Hz, which changes with the restoring forces 
of the spring couplings to a;/27r « 2.34 Hz (see App. 0. 

We harmonically excite one lattice site with a well- 
defined polarization by forcing the position of the two 
pendula. By tracking the position of all pendula one can 
extract the two-dimensional traces [xr,s(t),yr^s{t)]- From 


defines the wave vector in units of the inverse lattice con¬ 
stant. We show the resulting edge dispersions oj{k) in 
Fig. I^g). For each polarization, there is a chiral (that is, 
uni-directional) mode per gap, as expected. 

The edge dispersion indicates that the classical system 
faithfully implements the QSHE. However, the edge state 
of the QSHE are protected by the symmetry T. In the 
mechanical case, this symmetry is not related to time re¬ 
versal. Instead V is only symmetric under T if the ratios 
of the local couplings are correctly implemented. Hence, 
disorder in the spring constants induces both T-breaking 
terms that can remove the edge modes as well as benign 
disorder effects which only remove left- and right circular 
polarization as a global integral of motion. 

For a clean system, a linearly polarized excitation on a 
boundary site is split into two counter-propagating circu¬ 
larly polarized modes. Therefore the edge states can be 
used as a polarizing beam splitter. T-symmetry breaking 
diminishes the efficiency of the beam splitter by coupling 
the two polarizations. T-symmetric disorder, on the other 
hand, only leads to a selection of different polarizations in 
the splitting process. 

We measured the ratio of excitation ^r,s/ J2RAr,s 
between the left and the right long edge after exciting in 
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FIG. 4: Topological protection, a-d, Steady states with a 
sequence of removed sites illustrating the stability of the edge 
states against boundary roughness. The bottom right site is 
driven with linear polarization at a frequency in the upper 
gap. e, Wave packets launched at the marked site with lin¬ 
ear polarization at a frequency in the lower gap. Each clock 
represents an edge site on which the blue (red) wedge is cen¬ 
tered at the time to when the wave packet traverses the site. 
The angle of the wedge indicates the width a of the wave packet 
as illustrated in panel f. The color represents its polarization. 
From these clocks one can read off the propagation of the wave- 
packets throughout the system, f, Analysis of the wave-packet 
at a given site from which the passing time to and the width 
of the wave-packet a are extracted. 


the middle of a short edge (Fig.j^a)). We scanned all pos¬ 
sible polarizations by changing the relative amplitude and 
time-lag between the two local pendula. In analogy to the 
quantum mechanical case, we plot the splitting ratio as a 
colormap on a Bloch sphere, where the north pole (south 
pole) corresponds to right (left) circular polarization of the 
drive, respectively. The maximal imbalance between the 
left and right edge are not reached on the north (south) 
pole, but on two approximately antipodal points rotated 
by ~ 15° away from the poles (see Fig. [^b)). However, 


as illustrated in Fig. ic), we reach a splitting fidelity of 
99.80(4)% at the optimal points. From this we conclude 
that on the length-scale of our system, only disorder sym¬ 
metric under time reversal is relevant (see App. 11- 

Surface states in mechanical systems with co ncave 
boundaries are known as whispering gallery modes.^^^lHJ 
To demonstrate that our edge states are not mere whis¬ 
pering gallery modes and to highlight the robustness of the 
phononic edge states described here, we remove a sequence 
of sites from the dynamical problem, effectively creating 
a convex boundary. In Fig. |^a)-(d), we show the result¬ 
ing mode structure for a frequency in the lower gap. The 
result demonstrates that the exact shape of the boundary 
has no influence on the stability of the edge states, in con¬ 
trast to whispering-gallery modes, which are exquisitely 
sensitive to imperfections. 

To further strengthen the point that the edge states are 
topological rather than imposed by the finite-size geome¬ 
try, we create a domain wall between two different topo¬ 
logical sectors. We invert the effective flux seen by the 
two polarizations on six rows of the system (see App. 0. 
At the boundary between the two sectors, the spin Chern 
numbers change their values, which requires the presence 
of in-gap modes. We illustrate these topological in-gap 
modes along the sector boundaries by exciting a linearly 
polarized wave-packet on the short edge of the larger sec¬ 
tor. Thereby we also demonstrate the unique ability, ow¬ 
ing to the long timescales in our experiments, to inves¬ 
tigate dynamical rather than steady-state properties. As 
shown in Fig. Qe), the wave packet splits into two circu¬ 
larly polarized packets, each traveling along one edge. At 
the boundary to the second sector, they are deflected into 
the interior of the lattice (where they interfere to a linearly 
polarized wave-packet again), before they each travel in¬ 
dependently along the physical edge of the second sector 
in reversed directions. 

In this work, we have realized a mechanical ‘topologi¬ 
cal insulator’ displaying the phenomenology of the QSHE. 
Our results show that the topologically protected edge 
states are observable in a mechanical system despite the 
presence of disorder and decay processes. For poten¬ 
tial technological applications, the ability to define ‘edge’ 
states along domain-walls between different topological 
sectors is particularly promising. By stacking different 
domains, where the length of the domain walls can be con¬ 
trolled by the precise shape of the interface, stable acoustic 
delay lines can be engineered.Moreover, the work pre¬ 
sented here establishes the bridge between the quantum 
mechanical phenomena of topological insulators and me¬ 
chanical systems. The outlook to implement the acoustic 
counter-part of three-dimensional topological insulators is 
peculiarly interesting. 
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work. We acknowledge financial support from the Swiss 
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Appendix A: Dynamical matrix in real space 


Starting from Eq. ([^ we apply the unitary transforma¬ 
tion U to obtain D = U^PLU. The resulting equations of 
motion read explicitly 

^r,3s — (^0 T /*(^r±l,3s ^3r,3s±l); (-^1) 

yr,3s — (^0 ^f )yr,3s “t” f {jJrT.l,3s “t” y3r,3sTl')j (-^2) 

^r,3s-\-l — (^0 ^/*)^r,3s-t-l T /*^r,3s-t-l±l (-^^) 

/, ^ ^3/, ^ 

~ ■^(a;r'±l,3s-|-l) H-^(yr-|-l,3s-|-l — J/r-l.Ss-l-l); 


yr,3a+l — —(<^0 + ■^f)yr,3s+l + fyr,3s+l±l 


(A4) 


f, ^ V^// N 

2 (2/r±l,3s+lJ 2 ^r+l,3s+l ^r—l,3s+l); 


^r,3s+2 — (^0 -^^)^r,3s+2 “t“ f 


(A5) 


/ / N / X 

~ ^(2;r±l,3s-|-2) H-^( —2/r-|-l,3s+2 + 2/r-1.3s+2j) 


yr,3a+2 — ~{^0 + ^f)yr,3a+2 + f yr,3a+2±l 


(A6) 


f, ^ v^// ^ 

~ -^{yr±l,3a+2) H-^ (a^r-|-1.3s-|-2 — 1.3s+2)- 

Here, the term proportional to A = 3-|--\/3 encodes the 
restoring acceleration arising from the springs described 
by / = M/J, where J is the moment of inertia and M 
the torque constant of the springs. In our implementa¬ 
tion with pendula of 500 mm length and a mass of ap¬ 
proximately 500 g gravity gives rise to a bare frequency 
of a;o/27r « 0.75 Hz. The coupling constants are chosen 
to give \ff « 1.02 Hz. The above equations of mo¬ 
tions are implemented as outlined in Fig. For couplings 
with a positive sign in Eqns. (Al - A6[) in the s-direction 


simple springs between nearest neighbors are needed. For 
positive couplings in r-direction, two-lever arms connect¬ 
ing the pendula with springs achieve the same effect: A 
deflection of one pendulum induces a deflection of the sec¬ 
ond in the same direction. For negative couplings, this 
direction can be inverted via the use of only one lever arm 
as shown in Fig. We chose all lever arms and springs 
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light enough such that their inertia is negligible. The lat¬ 
tice spacing is 135 mm in r-direction and 120 mm in s- 
direction. Note that the effect of interchanging the flux 
$—>■—$ in "H amounts to an interchange of the sign in 
the cross-coupling between x and y pendula. Therefore, 
to create two different topological sectors one inverts these 
couplings on parts of the system. 


Appendix B: Band structure 


T-L is symmetric under time reversal T = iaylC, where ay 
is the Pauli matrix acting on the a-index and K, denotes 
complex conjugation. = — 1 which puts TL in symme¬ 
try class AII.ESl In two dimensions this class contains the 
QSHE and is characterized by a Z 2 topological index.l^ 
Both pseudo-spin sectors of H can be diagonalized and 
yield each three bands with dispersion 


)' (^r; 


2g27rz7/3 ^^^2/3 

/3{kr,ksy/^ 


where 


f3{kr,ks) = -cos{qr)-cos{qs) + \/[cos{qr) + COS {qs)y - 8 . 

Here 7 = 0,1,2 and qr{s) are the momenta in r (s) direc¬ 
tion, respectively. For our parameters the eigenfrequencies 
are given by u>-y{qk,qs) = \/s-y{qr: Qs) which are bounded 
by min[a;-y(gr-, Qs)] ~ 1.64Hz and max[w.y(( 7 r-, 9 s)] ~ 

2.91Hz. 


Appendix C: Spin Chern numbers 

In the pure case H describes two copies of the Hofs- 
tadter problem of a particle hopping on a square lattice in 
the presence of a flux $ = ^ per plaquette. Owing to the 
absence of any coupling between the two copies, a topo¬ 
logical index can be defined via the spin Chern numbers 
C = (C+jC-) evaluated for each sector separately.!^ The 
standarcP Z 2 index v is then given by 

ly = ^ mod 2. (Cl) 

In the present case the two copies a = ± have opposite flux 
and accordingly C+ = —C-. The Chern numbers for each 
gap can be calculated from the Diophantine equatiorP^ 
and yield C+ = I for the lower and C+ = — 1 for the 
upper gap yielding a non-trivial :/ = I. Any change in the 
Chern number is accompanied by a gap closing. Hence, the 
boundary between the different sectors of flux 4) = ±^, 
where C± change by ±2 carries two co-propagating edge 
channels per “spin”. 


Appendix D: Data analysis 


Time traces of all 270 pendula are measured in parallel 
with two cameras at an acquisition rate of 60 Hz. Due 
to the finite resolution of the cameras, the resulting posi¬ 
tion has an intrinsic error of 0.2 mm. Typical deflections 
are in the range of a few millimeter. This amounts to a 
relative measurement error between 10“^ and 10“^. More¬ 
over, the relation of the distance measured on the camera 
and the angle of deflection are in linear relation to a good 
approximation; we use mm throughout the text. We fit 
the obtained traces to the function 


f ^r,s(0\ 
\yr,s{t) J 


- At. 



A 


sin 

2 



/ sin(wt) 
\cos{u}t) 


+ cos 


1 

2 



I' Sin(-Wt -I- '&r,s)\ 1 
\^cos(-a;t -I- -dr,s)) j ' 


(DI) 


The resulting mean displacements A^ the polarizations 
ipr,s, and the angles 'dr,s have typically relative errors of 
maximally 10 “^ owing to the large number of measure¬ 
ments. Positive (negative) values of ipr,s correspond to left 
and right circular or elliptic polarizations. For ip^^s = 0; 
the angle 'dr,s corresponds the orientation of the linear po¬ 
larization. In the general case ipr,s and 'dr,s parameterize 
the polar and azimuthal angle of the Bloch sphere shown 
in Fig. § 


Appendix E: Localization length 


In order to assess the effect of T-symmetry breaking 
terms we calculated the localization length We used the 
KwantP^ package to determine the longitudinal conductiv¬ 
ity a{L) of H on a finite system of length L under the 
influence of disorder. The springs used in the experiments 
have an uncertainty of ^ 4% in their stiffness described 


by / in Fqns. (Al 


A 6 ). We assume Fqns. (Al A 6 ) to 


be disordered accordingly, transform back to the quantum 
mechanical problem H and extract ^ from the scaling be¬ 
havior of a{L) with system size. The resulting localization 
length depends on energy. For states in the gap we find 
the localization length of the edge states to be ^ « 1000 
lattice constants. For our system with a maximal distance 
along the edge of 46 sites, we therefore expect T-symmetry 
breaking to play no essential role for the edge channels. 
















